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Possible  Measurement  Errors 
In  Calibrated  AVHRR  Data 


1.  INTRODUCTION 

The  NOAA  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  is  a  scanning 
radiometer  that  simultaneously  senses  reflected  sunlight  and  and  emitted  thermal 
infrared  energy  in  the  five  “channels,”  or  window  regions,  listed  in  Table  1  (some 
AVHRR  instruments  do  not  contain  Channel  5;  NOAA-IO  is  the  last  AVHRR 
instrument  to  be  flown  without  Channel  5).  Channels  1  and  2  sense  reflected  solar 
energy  and  are  used  to  detect  cloud  cover,  snow  cover,  and  sea  ice,  along  with 
cyclones  and  volcanic  dust  plumes.  Clouds  reflect  sunlight  equally  well  in  both  the 
visible  Channel  1  and  the  near-infrared  Channel  2.  However,  land  backgrounds 
reflect  visible  light  poorly  when  compared  with  the  near-infrared.  Land  backgrounds 
that  appear  dark  in  the  visible  imagery  appear  bright  in  the  near-infrared  imagery. 
This  is  helpful  for  clear/cloud  distinctions  and  to  automated  cloud  analyses  such  as 
the  Air  Force  Global  Weather  Central’s  Real-time  Nephanalysis 
(AFGWC  RTNEPH)  because  most  backgrounds  will  have  a  different  spectral 
signature  than  clouds.  Thus  the  use  of  Channel  1  and  Channel  2  data  can  alleviate 
the  dependency  of  the  RTNEPH  on  background  brightness  fields. 

At  night,  the  infrared  Channel  3  sensor  of  the  AVHRR  measures  thermal  energy  in 
the  3.7p.m  wavelength  region.  Channel  3  was  chosen  to  complement  Channel  4  in  the 
remote  sensing  of  sea  surface  temperatures  by  providing  corrections  to  Channel  4 
temperature  measurements  contaminated  by  partial  cloud  cover  and  atmospheric 
water  vapor.  Atmospheric  water  vapor  is  a  weak  absorber  at  3.7p.m.  Consequently, 
thermal  energy  emitted  at  Channel  3  wavelengths  by  surfaces  such  as  land,  the  ocean, 


Table  1.  Spectral  Intervals  for  the  Five  NOAA-7  AVHRR  Window  Channels 


Five-Channel  NOAA-7  AVHRR 

Channel  1 
0.58-0.68p,m 

Channel  2 
0.725-1. lp.m 

Channel  3 
3.55-3. 93p.m 

Channel  4 
10.3-ll.3p.rn 

Channel  5 
ll.5-12.5p.rn 

and  clouds  can  penetrate  larger  amounts  of  water  vapor  without  being  attenuated. 

Thermal  energy  emitted  at  the  longer  Channel  4  wavelengths  can  not  penetrate  water 
vapor  as  effectively.  On  the  other  hand,  incident  solar  radiance  at  Channel  3 
wavelengths  is  not  negligible,  as  it  is  at  Channel  4  wavelengths.  Reflected  sunlight  is 
often  a  significant  part  of  daytime  Channel  3  radiance  measurements,  and  can  be  as 
large  as  emitted  terrestrial  radiances.  Channel  3  nighttime  imagery  has 
characteristics  similar  to  imagery  from  the  longer-wavelength  infrared  Channels  4 
and  5.  Channel  3  imagery  can  help  the  RTNEPH  to  discriminate  low  clouds  and  fog 
from  clear  land/ocean,  and  to  detect  thin  cirrus. 

Channel  4  senses  emitted  thermal  energy,  and  is  therefore  used  for  the  mapping  of  i 

thermal  radiation  emitted  by  clouds,  the  earth’s  surface,  and  oceans  during  both  day  I 

and  night.  Even  though  Channel  4  is  an  atmospheric  window  region,  radiation  j 

emitted  at  Channel  4  wavelengths  is  noticeably  affected  by  water  vapor.  Water  ( 

vapor  is  the  main  tropospheric  absorber  of  radiation  at  these  wavelengths,  and 
significantly  attenuates  emitted  thermal  radiation.  Channel  5  also  senses  emitted  i 

thermal  energy  and  has  the  same  characteristics  as  Channel  4,  except  that  Channel  5 
radiation  is  even  more  sensitive  to  water  vapor.  Corresponding  measured  brightness  j 

temperatures  from  both  channels  are  similar  for  clear,  cold  atmospheres  since  usually  1 

such  atmospheres  are  quite  dry.  But  due  to  attenuation  by  atmospheric  water  vapor, 

Channel  5  brightness  temperatures  are  as  much  as  3  percent  lower  than  their  , 

corresponding  Channel  4  temperatures  when  measured  in  moist  atmospheres. 

Channels  4  and  5  data  used  together  can  provide  cloud  analysis  models  such  as  the 

RTNEPH  with  atmospheric  water  vapor  attenuation  corrections  to  satellite-measured  I 

brightness  temperatures.  I 

It  is  the  intent  of  this  report  to  provide  AFGWC  with  information  on  how  AVHRR  ; 

data  are  measured  and  calibrated,  how  this  data  can  be  most  effectively  and  precisely 
used,  and  to  help  develop  an  improved  understanding  of  the  accuracy  and  precision  I 

that  can  be  expected  of  the  calibrated  data.  The  AVHRR  sends  down  10-bit  “raw  , 

count”  values  for  each  of  the  five  channels.  The  raw  visible  counts  (Channels  1  and  2) 
may  be  converted  to  normalized  albedoes,  and  the  raw  infrared  counts  (Channels  3 
and  4,  and  5  when  present)  may  be  converted  to  brightness  temperatures  using  , 

“calibration  coefficients.”  Calibration  coefficients  are  provided  for  each  channel  , 

individually  as  a  regular  part  of  the  data  stream,  and  can  change  every  Global  Area  I 

Coverage  scan  (GAC  scans  contain  4  km  data)  or  every  five  Local  Area  1 


Coverage/High  Resolution  Picture  Transmission  scans  (LAC/HRPT  scans  contain  1 
km  data).  Scan-to-scan  changes  are  most  noticeable  for  Channels  3,  4,  and  5.  A  brief 
description  of  the  procedure  for  converting  the  raw  counts  to  albedoes/temperatures 
is  given  in  Section  1.  The  variability  is  discussed  of  the  calibration  coefficients  on  a 
line-by-line  and  orbit-by-orbit  basis,  along  with  the  associated  effects  on  the  “look-up 
tables”  used  to  convert  from  counts  to  albedoes/temperatures.  In  Section  2  the 
consequences  of  truncating  10-bit  raw  counts  to  8-bit,  7-bit,  or  6-bit  counts  is 
discussed.  Section  3  summarizes  the  conclusions  of  Sections  1  and  2  and  discusses  the 
consequences  of  using  “constant”  look-up  tables  and  of  truncating  the  original  10-bit 
raw  data  stream. 


1 
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2.  AVHRR  CALIBRATION  PROCEDURES 


2.1  Visible  Data  Calibration 

Calibration  coefficients  are  provided  for  Channels  1  and  2  in  the  form  of  “slope” 
and  “intercept”  values.  Normalized  albedoes  are  a  linear  function  of  the  measured 
10-bit  count  value  and  of  these  slope  and  intercept  values.  However,  the  slope  and 
intercept  values  provided  must  first  be  scaled  by  230  and  222,  respectively,  before  they 
can  be  properly  used  to  convert  from  count  to  albedo. 

Having  scaled  the  slope  and  intercept  values,  they  may  now  be  used  to  convert  the 
AVHRR  visible  count  data  values  to  normalized  albedo.  The  percent  albedo  A  is 
computed  as  a  linear  function  of  the  count  value,  given  by 

Aj  n  mjC+bj,  (j  =  1,2),  (1) 

where  Aj  is  the  percent  albedo  measured  by  Channel  j;  m.j  is  the  scaled  slope  value  for 
Channel  j  in  units  of  percent  albedo/count;  C  is  the  measured  visible  10-bit  count 
value  for  Channel  j;  and  bj  is  the  intercept  value  for  Channel  j  in  units  of  percent 
albedo. 

Figure  1  contains  sample  plots  of  look-up  tables  for  converting  counts  to 
normalized  albedoes  for  Channels  1  and  2.  Note  that  as  the  count  C  increases,  the 
albedo  increases  as  well. 

For  the  NOAA-K,  -L,  and  -M  spacecraft  (scheduled  for  the  early  1990’s)  the 
normalized  albedo  look-up  table  will  be  piecewise  linear,  consisting  of  two  functions 
similar  to  Eq.  (1)  for  each  visible  channel.  One  linear  relationship  will  be  for  the 
lower  albedoes,  and  the  second  relationship  for  the  higher  albedoes;  they  will  be 
chosen  so  as  to  provide  the  finest  albedo  resolution  for  the  darkest  scenes. 
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2.2  Infrared  Data  Calibration 

Calibration  coefficients  are  provided  for  Channels  3  and  4  (and  5  when  present)  in 
the  form  of  “slope”  and  “intercept”  values.  Emitted  thermal  energy  (in  radiance 
units  of  mW/(m2  ster  cm-1))  is  a  linear  function  of  the  measured  10-bit  count  value 
and  of  these  slope  and  intercept  values.  However,  the  slope  and  intercept  values 
provided  with  each  scan  line  must  first  be  scaled  by  230  and  222,  respectively,  before 
they  can  be  properly  used  to  convert  from  count  to  radiance  and  finally  to  brightness 
temperature. 

Having  scaled  the  slope  and  intercept  values,  they  may  now  be  used  to  convert  the 
AVHRR  thermal  count  data  values  to  radiance.  The  radiance  is  computed  as  a  linear 
function  of  the  count  value,  given  by 


Ij  =  mjC+bj,  (j  -  3,4,5),  (2) 

where  Ij  is  the  emitted  upwelling  thermal  radiance  measured  by  Channel  j;  mj  is  the 
scaled  slope  value  for  Channel  j  in  units  of  mW/(m2  ster  cm-1)  per  count;  C  is  the 
measured  thermal  10-bit  count  value  for  Channel  j;  and  bj  is  the  intercept  value  for 
Channel  j  in  units  of  mW/(m2  ster  cm-1). 

Once  the  radiance  is  computed  for  thermal  Channel  j,  it  may  then  be  converted  to 
brightness  temperature  by  inverting  Planck’s  function.  The  Planck  function  relates 
monochromatic  radiance  with  wavenumber  (wavelength)  and  the  temperature  of  the 
emitting  surface.  It  can  be  written  as 

Cjl/3 

- ■  W 

eT  -1 

where  B  is  the  blackbody  Planck  radiance;  v  is  the  wavenumber  in  cm-1;  T  is  the 
temperature  of  the  emitting  surface  in  °K;  and  c2  and  c2  are  constants: 
Ci  =  1.1910659  x  10-6  mW/(m-2  ster  cm-4)  and  c2  =  1.438833  cm  °K.  Blackbody 
radiance  increases  with  temperature,  and  the  wavenumber  of  maximum  radiance 
increases  with  increasing  temperature. 

Inverting  Planck’s  function  for  temperature  yields 


T[B(.,T)|  - 


c2v 


In 


1+ 


B(v,T)  J 


(4a) 


Note  that  Eq.  (4a)  gives  a  form  for  temperature  as  a  function  of  radiance  B  and 
wavenumber  v.  Thus,  given  a  radiance  measurement  Ij(C)  (Eq.  (2))  for  thermal 


Channel  j,  the  equivalent  blackbody  brightness  temperature  can  be  computed  using  Ij 
in  place  of  B  in  (4a): 


T[Ij(C)j  = 


C2Vj 


In 

L  '■;>*  ] 

Ij(C)  J 

(4b) 


The  “central  wavenumber”  i>j  is  a  function  of  the  AVHRR  thermal  Channel  number 
j,  and  varies  from  satellite  to  satellite.  There  is  some  error  incurred  when  using  a 
constant  vy  due  to  details  of  the  sensitivity  of  the  sensor  itself  to  upwelling  thermal 
energy.  The  use  of  a  constant,  single  v-}  to  represent  the  complete  channel  j 
bandwidth  is  not  exact;  Eq.  (4b)  is  exact  only  when  Vj  is  allowed  to  vary  with 
radiance.  This  variation  is  not  large,  however,  being  is  greatest  for  Channel  3. 

Nominally,  for  the  NOAA-7  AVHRR,  v3  —  2671  cm-1  for  Channel  3,  v4  —  927 
cm-1  for  Channel  4,  and  i>6  =  841  cm-1  for  Channel  5;  these  are  the  values  used  for 
this  study.  More  detailed  information  on  each  Channel’s  central  wavenumber  is 
always  available  from  (the  most  recent  version  of)  the  NOAA  Polar  Orbiter  User’s 
Guide  (Kidwell,  1985).  It  is  important  to  remember  that  the  temperatures  obtained 
using  (4b)  are  not  corrected  for  atmospheric  attenuation. 

Equation  (4b)  can  be  rewritten  in  the  following  form  which  shows  brightness 
temperature  T  as  a  function  of  the  10-bit  count  C  explicitly.  Replacing  Ij  with  (2)  in 
Eq.  (4b)  yields 


Tj(C)  = 


In 


C2Vj 


cifj3 

mjC+bj 


(4c) 


where  mj  and  bj  are  as  for  Eq.  (2). 

Figure  2  contains  sample  plots  of  look-up  tables  for  converting  counts  to 
brightness  temperatures  for  Channels  3,  4,  and  5.  Note  that  as  the  count  C  increases, 
the  brightness  temperature  decreases. 


2.3  Variability  of  Calibration  Coefficients 

The  calibration  coefficients  can  change  every  GAC  scan  line,  as  mentioned  in  the 
introductory  section;  this  implies  the  look-up  tables  can  change  from  scan  to  scan  as 
well.  A  simple  statistical  analysis  was  performed  to  obtain  an  estimate  of  the 
variability  of  these  look-up  tables  on  a  scan-by-scan  basis.  Calibration  coefficients 
for  53497  GAC  scan  lines  were  sampled  from  the  NOAA-7  AVHRR,  representing  ~8 
half-orbits  in  the  northern  hemisphere  for  11  June  1982.  These  coefficients  were  then 


used  to  generate  look-up  tables  for  each  channel  on  a  scan-by-scan  basis. 

For  the  visible  channels,  the  mean  and  standard  deviation  was  computed  for  the 
53497  albedo-vs-count  look-up  tables  (one  look-up  table  for  each  scan  line).  For  the 
thermal  channels,  the  mean  and  standard  deviation  was  computed  for  the  53497 
brightness-temperature— vs— count  look-up  tables  (one  look-up  table  for  each  scan 
line).  The  average  curves  were  then  tabulated  along  with  their  “standard-deviation— 
envelopes.”  This  analysis  was  done  in  order  to  obtain  an  estimate  of  the  magnitude 
of  error  incurred  by  using  a  constant  look-up  table  to  convert  raw  counts  to 
normalized  albedoes  or  brightness  temperatures.  Results  are  discussed  in  the 
following  sections. 

2.3.1  VISIBLE  CALIBRATION  COEFFICIENTS 

For  each  of  the  53497  GAC  scan  lines,  a  look-up  table  as  given  by  Eq.  (1)  was 
calculated.  Then  the  mean  and  standard  deviation  of  the  albedoes  as  a  function  of 
10-bit  count  C  were  computed.  The  mean  is  defined  as 

<A,(C)>  =  2  A,, JC),  (5a) 

1N  n=l 

and  the  standard  deviation  as 

»A,(C)  -  [Aj,n(C)  -  <  A,(C)>]2  J“,  (5b) 

where  N  =  53497  GAC  scan  lines;  n  is  the  GAC  scan  line  number 

(n  =  1,2,3,...,53497);  Ajn(C)  is  the  albedo  look-up  value  for  10-bit  count  C,  GAC 
scan  line  n;  and  j  is  the  Channel  number  1  or  2. 

Figure  1  contains  plots  of  the  normalized-albedo— vs— counts  average  look-up 
tables  <Aj(C)>  for  Channels  1  and  2.  Table  2  contains  the  means  (from  (5a))  as  a 
function  of  count  for  Channels  1  and  2;  the  standard  deviations  aA,  (as  defined  by 
(5b))  for  all  counts  are  equal  to  zero  out  to  three  places  to  the  right  of  the  decimal, 
indicating  that  no  appreciable  error  is  to  be  expected  in  using  constant  look-up  tables 
for  the  “visible”  Channels  1  and  2. 


2.3.2  INFRARED  CALIBRATION  COEFFICIENTS 

For  each  of  the  53497  GAC  scan  lines,  a  look-up  table  as  given  by  Eq.  (4c)  was 
calculated.  Then  the  mean  and  standard  deviation  of  the  brightness  temperatures  as 
a  function  of  10-bit  count  C  were  computed.  The  mean  is  defined  as 


<Tl(C»  =  Tf  I  TjJC), 

1N  n=l 


(6a) 


Table  2.  Average  Look-up  Table  for  NOAA-7  AVHRR  Channel  1.  The 
mean  albedo  <A(C)>  is  given  by  Eq.  (5a)  for  j  =  1,2 


NOAA-7  AVHRR  Channels  1  and  2 
Normalized  Albedo  Look-up  Table 


Mean 

Albedo  (%) 
<A(C)> 


Channel  1 
10-bit  Count 
C 


Channel  2 
10-bit  Count 
C 
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where  N  =  53497  GAC  scan  lines;  n  is  the  GAC  scan  line  number 

(n  =  1,2,3,. ..,53497);  Tjin(C)  is  the  temperature  look-up  value  for  10-bit  count  C, 
GAC  scan  line  n;  and  j  is  the  Channel  number  3,  4,  or  5. 

Figure  2  contains  plots  of  the  brightness-temperature— vs— counts  average  look-up 
tables  <Tj(C)>  for  Channels  3,  4,  and  5.  Table  3  lists  the  means  and  Table  4  lists 
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and  the  standard  deviation  as 


ffTj(C)  -  Jt  [Tj,„(C)-<T,(C)»!  J  , 
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the  standard  deviations  (from  (6))  as  a  function  of  count  for  Channels  3,  4,  and  5. 
The  standard  deviations  crTj  are  not  zero  for  the  thermal  channels  as  they  are  for  the 
visible  channels,  indicating  that  some  error  can  be  expected  from  scan  to  scan  in 
using  a  “constant”  brightness-temperature  look-up  table,  i.e.,  a  look-up  table  that 
does  not  change  from  scan  to  scan  as  the  calibration  coefficients  do.  However,  the 
“smallness”  of  the  standard  deviations  listed  in  Table  4  indicates  that  any  errors 
incurred  are  most  often  small. 

Assuming  that  the  distribution  of  brightness  temperatures  to  counts  is  normal  for 
any  given  10-bit  count  value  C,  then  the  1-a  envelope  captures  67  percent  of  the 
brightness  temperatures  that  would  be  computed  from  count  C  using  Eq.  (4).  The  2 -a 
envelope  captures  95  percent  of  the  brightness  temperatures,  and  the  3-a  envelope 
captures  99  percent.  Thus  the  ±  1-a,  2-cr,  and  3-cr  envelopes  bound  the  errors 
incurred  67  percent,  95  percent,  and  99  percent  of  the  time  using  a  constant  look-up 
table.  It  follows  that  errors  greater  than  1-a  K  can  be  expected  in  3  of  every  10 
scanlines;  errors  greater  than  2 -a  K  in  1  of  every  20  scanlines;  and  errors  greater  than 
3-a  K  in  1  of  every  100  scanlines.  Errors  of  6-cr  K  occur  from  scan  to  scan 
Vioo  x  Vioo  (corresponding  to  a  deviation  of  3-(X  K  above  the  mean  in  one  scan  to 
3-a  K  below  the  mean  in  the  next  scan),  or  1  in  every  10,000  scanlines. 

Table  4  lists  the  magnitudes  and  corresponding  frequencies  of  errors  for  Channels 
3,  4,  and  5  that  can  be  expected  using  constant  brightness- temperature  look-up 
tables.  In  general  the  errors  for  Channel  3  will  be  the  largest  of  the  three  thermal 
channels  3,  4,  and  5;  this  is  related  to  the  fact  that  Channel  3  also  has  the  lowest 
signal-to-noise  ratio. 

Channel  3  1-a  errors  are  0.16  K  at  320  K,  decreasing  to  a  minimum  of  0.11  K  at 
287  K;  errors  then  increase  to  0.16  K  at  271  K,  0.20  K  at  267  K,  0.50  K  at  251  K,  and 
to  1.00  K  at  240  K.  After  240  K,  the  1-a  deviations  increase  rapidly  to  50  K  at  190 
K.  However,  at  3.7 p.m  the  emitted  190  K  radiance  is  • — 2  percent  of  the  emitted  320  K 
radiance.  For  this  reason  it  is  more  difficult  to  measure  3.7p,m  radiation  emitted 
from  very  cold  bodies  (since  the  signal  is  small  and  hence  the  signal-to-noise  ratio  is 
large),  which  in  turn  makes  it  unwise  to  rely  on  Channel  3  measurements  of  very  cold 
(190-230  K)  brightness  temperatures.  The  poor  accuracy  of  the  Channel  3 
measurements  for  very  cold  surfaces  causes  no  concern,  since  Channels  4  and  5 
measure  the  brightness  temperatures  of  such  surfaces  with  a  far  greater  certainty. 
Channel  3  is  most  accurate  in  measuring  higher  temperatures;  it  is  there  that  3.7|xm 
measurements  provide  the  most  useful  information  on  clouds. 

Both  Channels  4  and  5  have  much  lower  errors.  Channel  4  1-a  errors  decrease 
monotonically  from  0.12  K  at  320  K  to  0.06  at  190  K.  Channel  5  1-a  errors  are  even 
smaller,  decreasing  from  0.05  K  at  320  K  down  to  0.03  at  230  K,  and  then  rising  again 
to  0.05  K  at  190  K.  Thus  for  Channels  4  and  5  a  very  reasonable  estimate  of 
brightness  temperature  can  be  obtained  using  constant  look-up  tables. 
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NOAA-7  AVHRR  Channels  3,  4,  and  5 
Look-up  Table  Standard  Deviations 


Errors  (Standard  Deviations),  °K 


Channel  3 


Channel  4 

1-cr 

2-a 

3-a 

Channel  5 


Table  4.  Statistics  for  the  Look-up  Tables  of  the  NOAA-7  AVHRR  Thermal 
Channels  3,  4,  and  5,  Computed  Using  Eq.  (6b).  Errors  greater  than  ±  1-cr  °K  can  be 
expected  in  3  of  every  10  scaniines;  errors  greater  than  ±  2-a  °K  in  1  of  every  20 
scanlines;  and  errors  greater  than  ±  3-a  °K  in  1  of  every  100  scanlines 
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3.  EFFECTS  OF  TRUNCATING  THE  10-BIT  RAW  COUNTS  TO  8-,  7-,  AND 
6-BIT  COUNTS  (GRANULARITY) 

“Granularity”  is  defined  as  the  change  in  albedo  per  10-bit  count  for  the  visible 
data,  and  as  the  change  in  brightness  temperature  per  10-bit  count  for  the  infrared 
data.  The  granularity  in  albedo  or  in  brightness  temperature  can  be  kept  small  with 
10-bit  raw  counts,  since  10-bit  counts  allow  for  1024  possible  values  of  normalized 
albedoes  Aj  (which  range  from  0-100  percent)  or  1024  values  of  brightness 
temperatures  Tj  (which  range  from  190-320  K).  Granularity  increases  when  the 
number  of  bits  retained  in  the  raw  counts  is  decreased  (for  example,  in  going  from 
10-bit  down  to  8-,  7-,  or  6-bit  counts).  A  discussion  follows  of  the  effects  of  truncating 
the  AVHRR  10-bit  data  stream. 

3.1  Visible  Data  Granularity 

Mathematically,  granularity  for  AVHRR  visible  data  is  defined  as  the  first 
derivative  of  albedo  with  respect  to  10-bit  count.  Differentiating  Eq.  (1)  with  respect 
to  C  yields 

dA, 

<VlSL  =  n*  W 

where  GAj  is  the  granularity  in  units  of  percent-albedo  per  10-bit  count;  C  is  the  10- 
bit  visible  count;  and  rrij  is  the  scaled  slope  value  for  Channel  j  in  units  of  percent 
albedo  per  10-bit  count.  For  Channel  1,  mt  =  0.1068  percent-albedo/count,  so  that  a 
change  in  one  10-bit  count  corresponds  to  a  change  of  ~  0.1  percent-albedo.  For 
Channel  2,  m2  =  0.1069  percent-albedo/count,  so  that  a  change  in  one  10-bit  count 
likewise  corresponds  to  a  change  of  ~  0.1  percent-albedo. 

The  effect  of  reducing  the  raw  counts  from  10-bits  to  k-bits  is  simulated  by 
changing  the  scaled  slope  value  mj  as  follows: 

mj,k  =  210-k  mj,  (8a) 

where  mjk  is  the  scaled  slope  value  in  units  of  percent  albedo  per  k-bit  count;  k  is  the 
number  of  bits  retained  in  the  raw  counts;  and  mj  is  the  scaled  slope  value  for  10-bit 
counts,  for  Channel  j.  Replacing  Eq.  (8a)  into  Eqs.  (1)  and  (7)  gives 

Aj)k  =  mj|kCk  +  bj,  (j  =  1,2),  (8b) 
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where  GAj  k  is  the  granularity  in  units  of  percent-albedo  per  k-bit  count;  Aj  k  is  the 
percent  albedo  measured  by  channel  j;  mj(lc  is  given  by  Eq.  (8a);  Ck  is  the  k-bit  count 
value  which  ranges  from  0  to  2k — 1;  and  bj  is  as  for  Eq.  (1).  Note  that  for  k  —  10, 
Eqs.  (8b)  and  (8c)  reduce  to  Eqs.  (1)  and  (7),  respectively. 

Since  the  slope  mj  is  positive,  the  granularity  as  given  by  Eq.  (8c)  is  also  positive. 
This  is  because  albedo  increases  with  count.  It  is  also  easy  to  see  from  Eqs.  (8a)  and 
(8c)  that  the  granularity  for  the  visible  Channels  1  and  2  increases  by  a  factor  of  2  for 
each  bit  dropped  from  the  original  10-bit  raw  count  value  C.  Thus  for  Channels  1 
and  2  the  granularity  increases  from  0.1  percent-albedo  per  10-bit  count  to  0.4,  0.8, 
and  1.6  percent-albedo  per  8-bit  (k  =  8),  7-bit  (k  =  7),  and  6-bit  (k  =  6)  count, 
respectively. 

3.2  Infrared  Data  Granularity 

Mathematically,  granularity  for  AVHRR  thermal  data  is  defined  as  the  first 
derivative  of  brightness  temperature  with  respect  to  10-bit  count.  Differentiating  Eq. 
(4c)  with  respect  to  C  yields 


Gt  =  ' 

Tj  dC 


mjClc2vJ 


(1+-S?-)  I 

mjC+bj ' 


j  t  mjC+bj  J 


where  GTj  is  the  granularity  in  units  of  °K  per  10-bit  count;  C  is  the  10-bit  thermal 
count;  and  mj  is  the  scaled  slope  value  for  Channel  j  in  units  of  mW/(m2  ster  cm-1) 
per  10-bit  count.  All  other  variables  are  as  for  Eq.  (4c).  The  effect  of  reducing  the  raw 
counts  from  10-bits  to  k-bits  is  simulated  by  changing  the  scaled  slope  value  mj  as 
follows: 

mj,k  =  210-k  mj,  (10a) 

where  mj^  is  the  scaled  slope  value  in  units  of  mW/(m2  ster  cm-1)  per  k-bit  count;  k 
is  the  number  of  bits  retained  in  the  raw  counts;  and  mj  is  the  scaled  slope  value  for 
10-bit  counts,  for  Channel  j.  Replacing  Eq.  (10a)  into  Eq.  (4c)  gives 
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and  replacing  Eq.  (10a)  into  Eq.  (9)  gives 


dTjlk 

dCk 


mj,kclc2vj 


(10c) 


where  Gxjk  is  the  granularity  in  units  of  °K  per  k-bit  cox.  at;  Tj(k  is  the  brightness 
temperature  measured  by  channel  j;  mj)k  is  given  by  (10a);  Ck  is  the  k-bit  count  value 
which  ranges  from  0  to  2k— 1;  and  bj  is  as  for  Eq.  (2).  All  other  variables  are  as  for 
Eq.  (4c).  Note  that  for  k  =  10,  (10b)  and  (10c)  reduce  to  (4c)  and  (9),  respectively. 
Tables  5a-5c  list  the  granularities  for  Channels  3,  4,  and  5,  respectively,  as  a  function 
of  both  temperature  T  and  the  number  of  bits  k  of  the  raw  data  (that  is  10-bit,  8-bit, 
7-bit,  or  6-bit  counts).  Thermal  granularities  are  always  negative  because  mjk  is 
negative,  which  is  an  indication  that  brightness  temperature  decreases  with  increasing 
count.  However,  since  interest  lies  in  the  magnitude  of  the  granularities,  their 
absolute  values  are  listed  in  Tables  5-7  with  the  minus  sign  omitted  for  simplicity. 

As  can  be  seen  from  Eqs.  (9)  and  (10c),  and  from  Tables  5-7,  the  granularity  for 
the  thermal  channels  is  not  constant,  unlike  the  AVHRR  visible  channels.  The 
absolute  temperature  differences  increase  from  count  to  count  as  brightness 
temperatures  decrease.  This  is  because  AVHRR  thermal  sensors  are  designed  to  be 
most  accurate  in  measuring  radiances  emitted  by  warmer  surfaces.  The  exact  size  of 
the  differences  depends  on  the  count  C  (that  is,  on  the  temperature)  and  also  on  the 
AVHRR  channel  number  j.  For  example,  Table  5  shows  that  for  Channel  3  10-bit 
counts,  the  brightness  temperature  T3  changes  from  one  count  to  the  next  by  0.04  K 
at  310  K;  by  0.08  K  at  290  K;  by  0.18  K  at  270  K;  and  by  0.47  K  at  250  K.  Compare 
these  changes  with  those  in  Table  6  for  Channel  4  10-bit  counts  -  the  brightness 
temperature  T4  changes  from  one  count  to  the  next  by  0.07  K  at  310  K;  by  0.08  K  at 
290  K;  by  0.10  K  at  270  K;  and  by  0.13  K  at  250  K. 

Note  too  that  absolute  granularity  increases  as  bits  are  dropped  from  the  raw 
counts.  Although  the  granularity  is  not  constant  with  respect  to  counts  for  any  of  the 
thermal  Channels  3,  4,  and  5,  it  does  increase  by  a  constant  factor  of  2  for  each  bit 
dropped  from  the  original  10-bit  raw  count  value  C  (remember  this  occurs  for  the 
visible  channels  as  well).  This  is  seen  by  examining  (10c);  wherever  the  k-bit  count 
Ck  appears,  it  is  multiplied  by  the  weighted  slope  value  mj)k.  Remember,  Ck  has  been 
reduced  by  a  factor  of  210_k,  but  mjk  has  been  increased  by  a  factor  of  210-k.  Thus 
their  product  mj(kCk  remains  constant,  with  no  net  change  to  G-pj  k  wherever  this 
product  occurs.  However  in  the  numerator  of  (10c),  mjjk  appears  alone;  since  mjtk 
increases  by  a  factor  of  210-k  (see  (10a)),  it  causes  the  granularity  G-p,  k  to  increase  by 
a  factor  of  210-k  also.  For  example,  note  in  Table  7  that  for  Channel  5  the 
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Table  6.  Granularity  for  NOAA-7  AVHRR  Channel  4,  in  °K/Count,  as  a 
Function  of  Count  Sire  and  Temperature 


Temperature 

(°K) 

Channel  4  Granularity  (°K/Count) 

10-bits 

Count  Size 
8-bits  7-bits 

6-bits 

320 

0.06 

0.26 

0.52 

1.02 

310 

0.07 

0.28 

0.55 

1.11 

300 

0.08 

0.30 

0.61 

1.22 

290 

0.08 

0.33 

0.67 

1.35 

280 

0.09 

0.37 

0.74 

1.49 

270 

0.10 

0.41 

0.83 

1.66 

260 

0.12 

0.47 

0.93 

1.88 

250 

0.13 

0.53 

1.07 

2.14 

240 

0.15 

0.62 

1.23 

2.51 

230 

0.18 

0.73 

1.47 

2.93 

220 

0.22 

0.87 

1.74 

3.61 

210 

0.27 

1.06 

2.07 

4.36 

200 

0.33 

1.36 

2.63 

5.66 

190 

0.44 

1.79 

3.39 

6.79 

percent-albedo  per  count.  For  automated  nephanalysis  purposes  such  count-to-count 
albedo  granularities  more  than  likely  do  not  present  any  difficulty  to  the  quality  of 
the  visible  data.  (For  more  information,  refer  to  Section  2.1.) 


4.2  Infrared  Data 


Look-up  tables  for  the  AVHRR  thermal  Channels  3,  4,  and  5  vary  more  than  the 
visible  look-up  tables,  as  is  seen  from  Tables  3  and  4.  Although  the  errors  incurred 
by  using  a  constant  look-up  table  are  usually  small,  they  can  become  appreciable  in 
size  from  time  to  time  as  previously  discussed;  these  errors  are  largest  for  Channel  3. 
Whether  or  not  these  errors  are  significant  depends  on  the  application  of  the  infrared 
data.  For  nephanalysis  purposes  the  temperature  deviation  errors  incurred  using 
constant  look-up  tables  for  the  longwave  Channels  4  and  5  will  most  likely  not 
present  any  difficulty  to  the  RTNEPH  in  terms  of  the  quality  of  that  infrared  data. 
However,  Channel  3  constant— look-up— table  errors  can  be  more  significant, 
especially  for  brightness  temperatures  270  K  and  lower;  it  is  for  this  range  of 
temperatures  that  the  Channel  3  sensor  is  most  inaccurate,  often  requiring  updated 
calibrations  on  a  scan-by-scan  basis  (See  Table  4).  In  comparison,  for  multispectral 
techniques  that  are  sensitive  to  interchannel  differences  in  brightness  temperatures, 
look-up— table  errors  can  be  significant  for  Channels  4  and  5  as  well.  For  example, 
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Table  7.  Granularity  for  NOAA-7  AVHRR  Channel  5,  in  °K/Count,  as  a 
Function  of  Count  Size  and  Temperature 


Channel  5  Granularity  (°K/Count) 

Temperature  Count  Size 

(°K)  10-bits  8-bits  7-bits  6-bits 


algorithms  using  multispectral  infrared  measurements  estimate  water  vapor  amounts 
using  the  ratio  of  Channel  4  and  Channel  5  brightness  temperatures;  this  ratio  can  be 
quite  sensitive  to  measurement  errors  of  tenths  and  even  hundredths  of  degrees  K. 
(For  further  description  refer  to  Sections  1.2  and  1.3.2,  Tables  3  and  4,  and  Figure  2.) 

If  the  thermal  data  is  truncated  from  10  bits  to  8,  7,  or  6  bits  then  several 
consequences  arise.  First,  the  granularity  throughout  the  temperature  scale  will  be 
degraded.  In  going  from  10-bit  to  6-bit  counts  (6-bit  counts  are  what  the  AFGWC 
Satellite  Global  Data  Base  (SGDB)  currently  retains),  the  brightness  temperature 
granularity  increases  by  a  factor  of  16.  But  more  importantly,  the  accuracy  of  the 
brightness-temperature  differences  among  the  three  thermal  channels  can  also 
decrease  by  a  factor  of  16  or  more.  The  thermal  data  may  be  so  degraded  that  it 
would  be  difficult  for  multispectral  algorithms  to  use  it  effectively.  The  success  of 
such  algorithms  is  often  dependent  on  small  differences  in  brightness  temperature 
from  one  thermal  channel  to  the  next.  (Refer  to  Section  2.2  and  Tables  5a-5c.) 
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4.3  General  Comments 

Finally,  note  that  the  temperature  look-up  tables  (see  Table  3)  are  not  the  same 
among  the  AVHRR  thermal  Channels  3,  4,  and  5.  It  should  also  be  noted  that  the 


visible  Channels  1  and  2  albedo  look-up  tables,  although  nearly  the  same  for  NOAA- 
7,  need  not  necessarily  be  so  similar  for  other  NOAA  polar  orbiters.  In  other  words,  a 
visible  raw  count  may  represent  different  normalized  albedoes.  Similarly  a  given 
thermal  raw  count  value  represents  quite  different  temperatures  depending  on  the 
channel  in  which  it  was  measured. 

For  this  reason,  if  it  is  desired  to  intercompare  or  combine  data  from  two  or  more 
of  these  channels  into  a  single  image  product  a  two-step  process  must  take  place. 
First,  the  raw  counts  must  be  converted  to  brightness  temperatures  (or  albedoes,  if 
data  from  the  visible  channels  is  to  be  combined)  using  the  appropriate  look-up 
tables.  Next,  these  brightness  temperatures  (albedoes)  can  be  converted  to  a  common 
grayshade  scale  for  image  processing  purposes.  Always  remember  that  it  is  incorrect 
to  combine  uncalibrated  “raw  count”  data  from  two  or  more  of  the  thermal  channels 
(or  from  the  two  visible  channels)  without  first  converting  the  raw  counts  to 
brightness  temperatures  (or  normalized  albedoes). 

In  the  final  analysis,  the  errors  incurred  using  constant  look-up  tables  or  in 
truncating  the  original  10-bit  raw  count  data  stream  will  be  tolerable  for  some 
applications  and  not  for  others.  The  extent  to  which  a  high  accuracy  is  needed  for 
measurements  of  normalized  albedoes  and  brightness  temperatures  wil’  be  driven  by 
individual  algorithms  and  applications.  In  general,  future  multispectral  techniques 
that  use  thermal  data  will  most  likely  require  higher  temperature  measurement 
accuracies  than  the  current  6— bit  SGDB  can  support.  It  is  strongly  recommended 
that  the  AFGWC  retain  at  least  7  and  preferably  8  bits  of  thermal  data  in  both  the 
Satellite  Data  Handling  System  (SDHS)  and  the  SGDB. 


Appendix  A 

Average  NOAA-7  Slope  and  Intercept  Calibration  Coefficients 


The  average  NOAA-7  AVHRR  calibration  coefficients  for  the  53497-line  sample 
used  in  this  study  are  tabulated  below: 

AVHRR  Calibration  Coefficients  For  10-bit  Raw  Count  Data 


Channel 

Slope  m 

Intercept  6 

1 

0.106806 

-3.4399 

2 

0.106905 

-3.4882 

3 

-0.001477 

1.4680 

4 

-0.154670 

152.2701 

5 

-0.180880 

177.9360 

For  the  visible  Channels  1  and  2,  the  listed  slope  can  be  used  for  mj  in  Eqs.  (1),  (7), 
and  (8a);  the  listed  intercept  can  be  used  for  bj  in  Eqs.  (1)  and  (8b).  The  visible  slope 
values  are  in  units  of  percent-albedo  per  10-bit  count.  The  visible  intercept  values  are 
in  units  of  percent-albedo. 

For  the  thermal  Channels  3,  4,  and  5,  the  listed  slope  can  be  used  for  mj  in  Eqs. 
(2),  (4c),  (9),  and  (10a);  the  listed  intercept  is  bj  as  for  Eqs.  (2),  (4c),  (9),  (10b),  and 
(10c).  The  thermal  slope  values  are  in  units  of  mW/(m2  ster  cm-1)  per  10-bit  count. 
The  thermal  intercept  values  are  in  units  of  mW/(m2  ster  cm-1). 
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